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ORTHOGONAL  DECOMPOSITION  OF  FLOW  FIELDS  IN  DISTANCE/DISTANCE/TIME  SPACE 
INTRODUCTION 

Bethke  and  Viets  have  developed  two  techniquea  for  identification  of 

fl  21 

rotational  structures  (vortices)  in  flow  fields.  ’  The  first  utilizes 
phase  information  in  the  discrete  Fourier  transform  of  the  field  while  the 
second  involves  spatial  behavior  of  local  orthogonal  decomposition  components 
of  the  field.  While  the  second  method  is  generally  superior  to  the  first, 
both  methods  required  whole  field  input  data  at  a  given  instant  in  time.  This 
research  addresses  the  relaxation  of  this  restriction  through  the  analysis  of 
local  field  behavior  as  a  function  of  time,  and  is  based  on  the  second  method. 
ORTHOGONAL  DECOMPOSITION  IN  THE  DISTANCB/DISTANCB/TIMB  DOMAIN 

As  an  approach  to  distance/distance/time  (d/d/t)  orthogonal  decomposition, 
suppose  we  consider  flow  velocities  taken  at  4  fixed  rectangularly  related 
locations  at  two  different  times.  These  are  depicted  in  Figure  1  as 
velocities  1-4  at  one  sample  time  say  k-1,  and  velocities  5-8  at  the  next 
sample  time,  k.  The  d/d/t  data  cell  contains  8  "locations"  each  having  a  u 
and  v  spatial  velocity  component.  Thus  at  any  two  adjacent  sample  times  we 
have  a  data  cell  with  16  velocities,  Uj,v,,  through  u, ,  vf.  The  approach 
adopted  here  will  model  the  behavior  of  the  flow  in  a  spatially  fixed  data 
cell  and  then  study  the  cell  behavior  as  a  function  of  time  to  extract  flow 
field  information. 
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Figure  1.  Distance/Distance/Time  Flow  Data  Cell 


The  data  in  a  d/d/t  flow  cell  will  be  aodeled  with  12  independent  measures 
of  its  behavior.  These  12  descriptors  (a’s)  will  describe  the  local  time 
behavior  of  six  spatial  flow  measures  (fi’s)  at  the  two  different  cell  times  (k 
and  k-1).  A  set  of  six  fi'a  describe  the  amount  of  rotation,  horizontal  and 
vertical  translation,  horizontal  and  vertical  expansion,  and  shear  at  each 
sample  time.  That  is  at  a  given  sample  time,  say  k-1,  the  velocities 
ui » vi >u2 » v2 »ui » vi »u* * v«  wiH  be  described  by  a  0  weighted  set  of  basis  vectors 
plus  an  error  vector  (c’s),  i.e. 
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It  is  seen  that  the  0  basis  vectors  describe  the  six  flow  types  shown  in 
Figure  2,  where  A  and  B  are  the  vertical  and  horizontal  dimension  of  the  data 


cell. 
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ROTATION  TRANSLATION 


TRANSLATION  EXPANSION 


EXPANSION  SHEAR 


Figure  2  0  Basis  Vector  Components 

The  data  in  a  d/d/t  flow  cell  will  be  modeled  by  least  squares  fitting  the 
6  orthogonal  descriptors,  /J’s,  at  each  of  the  two  sample  times.  At  sample 
time  k-1  this  will  give  /J, (k-1) , . . . ,0f (k-1)  and  at  sample  time  k, 
0X  (k) , . . . ,0t (k) .  Time  will  be  incorporated  via  the  sum,  aia(k),  and 

difference,  a^(k) ,  of  the  0  types.  That  is, 

au  (k)  =  £[/>,<k)  *  />i(k-l>] 

and 

ai(J(k)  =  £[/>.(k)  -HjCk-l)]. 
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Thus  each  distance/distance/time  data  cell  will  be  described  by  the  12  a 


parameters  . . .  •ag8»ai<j*  •  •  •  »asd* 

The  modeling  process  described  above  can  be  viwed  as  a  linear 
squares  operation  in  which  the  vector  composed  of  the  16  velocities  of  a 
data  cell,  u, ,Vj , . . .ua , vt ,  is  to  be  modeled  as  an  a  weighted  combination 
basis  vectors  plus  an  vector  of  errors,  4,  That  is, 
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sample  time  k  sample  time  k-1 


or  in  matrix  notation 


V  =  X  a  +  t 

Note  that  the  1st  and  7th  basis  vectors  correspond  to  the  tine  sun  and 
difference  of  spatial  rotation  respectively,  the  2nd  and  8th  vectors  the  tine 
sun  and  difference  of  spatial  horizontal  translation,  and  so  forth.  Also  note 
that  the  basis  vectors  are  all  nutually  orthogonal  which  not  only  greatly 
sinplifies  least  squares  estimation  of  the  a  values  but  makes  these  estimates 
statistically  independent. 

/K 

Least  squares  estimates  of  the  a  values,  a,  are  readily  found  using 

S  ’  t*TX1_1  xTv  =  l  XTV 

2  2 

where  A  +B  =  1. 

PROPERTIES  OF  THE  TIME  DOMAIN  SUM  AND  DIFFERENCE 

As  noted  above,  this  research  is  based  on  tine  domain  behavior  of  a  given 

orthogonal  flow  discriptor,  (/T),  as  exemplified  by  the  sum  and  difference  of 

this  discription  at  adjacent  sanple  times.  That  is,  for  the  ith  orthogonal 

descriptor  0.  the  sum  descriptor,  a.  at  sample  times  k  and  k-1,  is 
i  is 

ais(k)  =  j(/5.(k)  +  /J.(k-l)) 

or  more  conviently  (i  understood) 

ag(k)  =  j(/J(k)  +  /J(k-l)). 

The  difference  descriptor  is  thus  written  as 

ad(k)  =  ^(/J(k)  -  /S(k-l) ) . 

These  sum  and  difference  operations  can  be  viewed  as  time  domain  digital 
filtering  of  the  0  vlaues  to  produce  the  a  values  as  depicted  in  Figure  3. 


Figure  3.  Tiae  domain  filtering  operations 


Insight  into  the  sua  and  difference  operations  can  be  gained  through  their 
frequency  response.  For  exaaple,  if  we  take  the  Z-transfora  of  the  sua 
operation 

as(X)  =  £(4(K)  +  /*(  K-1)) 

we  get 

as(Z)  =  2(4(Z)  + 

The  ratio  of  the  output  and  input  produces  the  transfer  function,  H  (Z),  of 
this  filtering  operation 


aCZ)  1  -i  7 

Vz)  ■  -5m  3  5a+z  1 3  - 

jwT 


(Z+l) 
In¬ 


setting  the  complex  variable  Z  equal  to  eJ  ,  where  T  is  the  tiae  between 

saaples  we  get  the  frequency  response  of  the  filtering  operation,  i.e. 

I,  jwT^.v 

„  ,  J«*T,  _  5(e  11 

“.<•  >  * - 3=T- 


1  ej»T/2jejwT/2  +  e~ j«T/2  j 


The  difference  operation  likewise  has  the  transfer  function 

i(z-l) 

Hd(z)  =  -r— 

and  the  frequency  res ponce 

Hd(e^T)  =  sin(«T/2) . 

The  frequency  response  Magnitudes  for  the  sum  and  difference  oeprationa 
are  shown  in  Figure  4. 


Figure  4.  Frequency  Response  Magnitudes 
It  is  seen  that  the  sum  and  difference  are  low  and  high  pass  filter  operations 
with  gains  of  cos(wT/2)  and  sin(wT/2)  respectively.  That  is,  the  sun 

operation  enhances  the  low  while  suppressing  the  high  frequency  tenporal 
behavior  of  a  given  flow  descriptor.  The  difference  operation  does  just 

the  opposite,  suppressing  low  and  enhancing  high  tenporal  frequencies.  It 


should  be  noted  that  "low”  and  "high"  frequencies  here  are  relative  to  one 


half  the  sample  frequency,  the  highest  unaliased  frequency  in  the  data,  ir/T 
radians/ second . 


The  sin  and  difference  operations  cure  also  related  to  the  discrete  Fourier 

Transform  (DFT)  which  transform  N  discrete  tiae  data  values  to  N  discrete 

frequency  (spectral)  values  via  the  operation 

N-l 

p  f  mr%\  -  \ 


F(«2)  =  J  f(KT)e" 


where 


For  the  two  point  case  (N=2)  this  becomes 


F(Ofi)  =  f(OT)  +  f (IT) 


F(1I2)  =  f(0T)  -  f(lT)  . 

which  is  proportional  to  the  previously  discussed  sua  and  difference 
operations  but  with  the  term  reversed  in  the  difference.  Repetative  use  of 
sua  and  difference  to  describe  temporal  behavior  of  a  given  flow  parameter, 
0^,  can  thus  be  considered  as  akin  to  repetitively  computing  its  2-point 
frequency  spectrua  as  a  function  of  tiae. 

SYWTHBTIC  VORTEX  TEST  RESULTS 

The  synthetic,  counter  clockwise  velocity  field  shown  in  Figure  5  was  used 
to  teat  the  orthogonal  decomposition  method.  The  field  was  assumed  to  be 
moving  to  the  right  at  a  constant  velocity  with  unchanging  rotational 
behavior.  The  distance/distance/tiae  data  cells  were  formed  using  the 
velocities  at  4  rectangularly  adjacent  locations  to  respreaent  a  sample  tiae, 
say  k-1,  and  the  4  location  velocities  imediately  to  the  left  as  tiae  k. 


Figure  J.  Velocity  distribution  in  the  synthetic 
vortex. 


Bach  of  the  12  a  parameters  was  computed  for  each  of  the  d/d/t  data  cells,  so 
constructed,  in  the  field. 


The  synthetic  field  results  are  displayed  in  a  set  of  6  vector  field 
plots,  one  plot  for  each  flow  descriptor  (0)  type;  rotation,  horizontal 
translation,  and  so  forth.  Each  vector  in  a  given  plot  is  formed  by  the  sum 
and  difference,  a  and  a.,  values  of  a  d/d/t  data  cell.  The  sum,  a  ,  and  the 

9  O  9 

difference,  a,,  being  the  horizontal  and  vertical  vector  component  lengths 
d 

respectively.  The  origin  of  each  vector  is  determined  vertically  by  the 
vertical  location  of  the  data  cell  in  the  flow  field  and  horizontally  by  the 
data  cell  sample  time,  k.  To  readily  allow  comparison  with  smoke  pictures  and 
other  instantaneous  time  field  deacriptions,  the  earliest  time  data  is  on  the 
right  of  the  plots,  with  time  increasing  to  the  left. 

These  plots  exhibit  the  behavior  one  would  see  monitoring  the  velocities 
at  4  rectangularly  oriented  locations  as  a  funciton  of  the  flow  descriptor 
type  (0  plot  type),  vertical  cell  location  (plot  row),  and  time  (location 
within  the  plot  row);  as  exhibited  by  the  a  vectors. 

The  synthetic  vortex  0t  plot  (rotation),  Figure  6,  shows  vortex  activity 
along  a  given  row  as  viewed  right  to  left  (increasing  time).  That  is,  the 
horizontal  length  of  the  vector  grows  and  subsides  indicating  vortex  strength, 
and  the  vertical  length  goes  positive  and  then  negative  indicating  increasing 
and  then  decreasing  vorticity  as  the  vortex  passes.  The  time  trajectory  of 
the  vector  along  a  given  row  thus  indicates  the  presence  and  rate  of  change  of 
the  vorticity  at  that  vertical  location  in  the  flow  field.  The  characteristic 
vector  time  trajectories  for  rows  above  and  below  the  vortex  center  are 
displayed  on  the  plot  for  this  counter  clockwise  example.  A  like  field, 
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He  hor'  ZoH*-*. 

rotating  clockwise,  would  caus elector  components  dle{flley«<fM*er«'  to  be  of 
opposite  sign  at  all  times  and  would  result  in  the  dashed  trajectories. 

The  4^  (rotation)  plot  clearly  indicates  the  presence  of  vorticity  that 
is,  having  vectors  of  significant  non-zero  length.  The  horizontal  center  of 
the  vortex  is  indicated  by  a  horizontal  vecotr,  i.e. ,  one  with  vorticity 
but  of  unchanging  strength.  The  vertical  center  of  the  vortex,  however,  is 
difficult  to  locate  using  only  one  vector  row  (set  of  4  sensors)  because  each 
row  has  a  similar  vector  trajectory.  Observation  of  the  ^(rotation)  vector 
thus  indicates  the  presence  of  vorticity  and  the  streaawise  location  of  the 
vortex  center  but  not  the  anti-streamwise  center  location. 

Figure  7,  the  horizontal  translation  plot,  provides  infornation  not 
present  in  the  Figure  6  rotation  plot.  The  vector  trajectoris  are  all 
clockwise  (c.w.)  but  the  sun  (horizontal)  vector  component  (a  )  are  negative 

3 

above  the  vortex  center  and  positive  below.  A  singlerow  in  Figures  6  and  7 
can  thus  locate  the  vortex  center  in  tine  and  indicate  if  the  rotational 
center  is  above  or  below  the  observation  row.  Figure  8,  the  vertical 
translation  plot,  again  shows  c.w.  vector  trajectories  but  in  this  case  the 
trajectories  go  both  positive  and  negative  in  the  horizontal  direction.  These 
trajectories  are  of  similar  nature  above  and  below  the  vortex  center,  and  are 
thus  of  use  in  detecting  vorticity,  its  rotational  direction,  and  locating  it 
in  tine,  but  not  in  the  antistreanwise  direction. 

Figures  9  and  10,  horizontal  and  vertical  expansion,  each  show  different 
patterns  above  and  below  the  center  streanline,  allowing  the  vortex  center  to 
be  so  identified  as  well  as  in  time.  Rotation  direction  is  also  evident. 

Figure  11,  shear,  produces  unique  double  rotation  trajectories  suitable 
for  identification  of  the  vortex  center  in  the  streanwise  but  not  the 


10 


ant is t 


ise  directions. 


The  c.w.  trajectories  indicate  c.c.w.  vortex 


rotation. 

D/D/T  DECOMPOSITION  OF  AW  ACTUAL  FLOW  FIBLD 

The  analysis  and  display  techniques  applied  to  the  synthetic  flow  field 
were  applied  to  the  field  shown  in  Figure  12.  This  field  was  generated  by  a 
can  shaped  rotor,  rotating  at  the  wall  of  a  low  speed  wind  tunnel. 


Figure  \J. 


Smoke  flow  visualization  of  the  rotor 
flowfield. 


A  set  of  instantaneous  flow  field  velocities  froa  the  tunnel  are  shown  in 
Figure  13. 


Figure  13  Instantaneous  flowfield  velocities 
The  instantaneous  velocity  field  shown  in  Figure  13  was  analyzed  as  time 
invariant  translational  data  using  the  same  techniques  as  above.  Figure  14 
displays  the  rotational  (/8^)  behavior.  The  clockwise  (c.w.)  vortices  in  the 
lower  left  and  upper  right  are  clearly  apparent  as  indicated  by  the  large 
horizontal  vector  components  pointing  to  the  left.  A  small  c.c.w  vortex  is 
also  noted  at  the  lower  center  of  the  field  as  indicated  by  short  vectors  to 
the  right.  It  is  seen  that  any  one  row  of  data,  4  senor  locations  as  a 
function  of  tine,  provides  vorticity  information  along  that  row  as  well  as  the 
streanwise  locations  of  the  vortex  centers.  In  addition  it  is  seen  that 
vector  behaviors  along  rows  in  or  near  vortex  centers  exhibit  strong  smooth 
behaviors  while  those  along  the  vortex  edge  exhibit  more  erratic,  high 
frequency  movements. 

12 
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Figure  14  shows  the  time  behavior  of  the  horizontal  translation  component 
of  the  field.  It  is  seen  that  all  vecotrs  in  the  field  have  a  positive  bias 
direction  indicating  general  field  flow  to  the  right.  Vortex  activity  along  a 
row  is  therefore  best  indicated  by  the  vertical  activity  of  these  vectors.  In 

the  rows  below  a  c.w.  vortex  center  the  vectors  nod  down  and  then  up  with 

time.  In  the  rows  above  they  go  up  and  then  down.  Given  the  fact  these 
vortices  are  c.w.  (from  Figure  180  this  information  not  only  indicates  the 
streamwise  location  of  the  vortices  but  also  if  they  are  above  or  below  the 

observed  row.  It  is  also  interesting  to  note  that  the  vectors  in  this  field 

£  u 

seem  to  generally  depict  the  trajectory  of  the  vortex  center.  Figure  if  shows 
the  time  activity  of  the  vertical  translation  activity  parameter  />, .  As 
expected,  vortex  activity  is  indicated  by  strong  sweeping  movements  along  a 
given  row.  The  streamwise  location  of  the  vortex  being  indicated  being 
vertical  down  vector  (c.c.w)  or  vertical  up  vector  (c.w.)  It  is  noted  that 
along  any  row  the  vectors  rotate  approximately  the  same  amount,  (one 
revolution  c.w.  independent  of  the  row). 

n  17 

Figures  and  show  the  cell  horizontal  and  vertical  expansion  time 

activity  along  a  row  (time,  as  usual,  increasing  right  to  left).  Comparison 
of  these  figures  with  their  synthetic  vortex  counterparts,  Figures  9  and  10, 
once  again  indicates  the  existence  of  and  streamwise  location  of  the  vortex, 
its  rotation  direction,  and  if  the  vortex  center  is  below  or  above  the 
observed  row. 

Figure  1£  shows  the  row  time  behavior  of  the  shear  parameter,  fit.  The 
characteristic  double  oscillation  behavior  shown  in  the  synthetic  case  is 
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again  seen.  As  before,  the  regions  of  vorticity  are  narked  by  significant 
activity  in  the  tine  behavior  of  this  vector. 

CONCLUSION 

This  study  explores  the  utility  of  locating  rotational  structures  in  a 
flow  field  by  locally  deconposing  the  fields  spatial  and  tine  behavior  via  a 
set  of  orthogonal  descriptors  and  then  viewing  these  descriptions  as  a 
function  of  tine.  It  is  seen  that  these  descriptions,  based  on  a  8  point 
distance/distance/tine  cell,  can  identify  rotational  structures,  their 
direction  of  rotation,  their  streanwise  location,  and  if  they  are  above  or 
below  the  streanwise  observation  location. 

The  behavior  of  the  cell  paraneters,  as  a  function  of  tine  (streanwise) 
and  antistreanwise  displacement,  was  very  distinct  and  clearly  exhibited  a 
great  deal  of  structure  in  addition  to  those  aspects  discussed  above.  It  is 
reasonable  that  the  nethod  of  flow  decomposition  and  display  presented  here 
nay  offer  additional  insights  into  observed  flow  fields  dynanics. 

Further  investigation  might  be  conducted  via  the  four  lcoation  two  sample 
tine  cell,  as  above,  or  perhaps  by  a  two  location  four  sanple  tine 
configuration.  In  either  case  the  tine  behavior  of  the  computed  vectors  night 
well  be  displayed,  recorded  and  studied  via  a  storage  osillosope. 
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Abstract 

Forced  unsteady  flows  are  examined  from  the 
point  of  view  of  potential  application  of  new 
devices  which  are  made  possible  by  the  improved 
understanding  of  chese  flows.  In  particular, 
flows  are  examined  which  require  no  external 
driver  but  obtain  their  energy  from  the  free- 
stream  velocity.  Other  flows  require  no  moving 
parts  at  all  to  generace  the  unsteadiness. 

Introduction 

Unsteady  separated  flows  may  appear  as  a  con¬ 
sequence  of  the  motion  of  the  surface  of  the  body 
or  as  a  result  of  the  unsteadiness  in  the  flow  over 
thac  surface.  In  many  cases,  when  confronted  by 
such  flows,  the  objective  is  to  eliminate  or  at 
least  control  the  separation. 

The  currenc  research  program  has  taken  a  some¬ 
what  different  perspective.  The  objective  of  the 
program  is  to  examine  the  dynamics  of  flows  which 
are  forced  to  become  unsteady  and  to  specifically 
consider  the  possibilities  for  the  application  of 
such  flows  as  well  as  determining  their  basic 
structure.  In  addition  to  the  question  of  identify¬ 
ing  devices  which  can  drive  a  flow  unsteady  in  a 
suitable  way,  it  is  necessary  to  determine  the  com¬ 
position  and  dynamics  of  the  large  scale  structures 
creaced  in  the  periodic  flow.  The  latter  point  has 
led  to  techniques  which  have  some  utility  in  deter¬ 
mining  che  vortex  structure  of  quasi-steady  three 
dimensional  flows. 

When  considering  possible  techniques  to  pro¬ 
duce  an  unsteady  periodic  flow  from  a  nominally 
steady  flow,  two  basic  requirements  are  apparent: 

a.  The  device  must  be  mechanically  simple 
and  easy  to  build; 

b.  [f  possible,  the  device  should  not  require 
external  power  to  drive  it. 

Past  Efforts1  7 


flow  visnalizatlon  results  of  Figure  1)  which  have 
an  important  effect  on  the  flowfield.  Although 
these  nozzles  can  be  employed  to  control  boundary 
layers,  most  of  the  effort  has  been  applied  to  free 
jets.  The  jet  is  controlled  by  the  fluidic  feed¬ 
back  system  wrapped  around  the  body  of  the  jet 
nozzle  in  Figure  2.  An  alternate  version  employing 
a  different  fluidic  feedback  has  been  found  to  be 
capable  of  producing  an  unsteady  jet  even  if  the 
phase  of  the  jet  flow  is  different  from  that  of  the 
surroundings. 


Figure  1.  Smoke  flow  visualization  of  a 
fluidically  oscillating  jet. 
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Figure  2.  Schematic  of  the  fluidically 
oscillacing  jet  nozzle. 


Under  che  present  research  program,  two 
devices  have  been  employed  to  generace  che  time 
dependence,  the  first  is  a  series  of  purely  fluid- 
ically  controlled  Jet  nozzles  which  satisfies  both 
of  the  requirements  cited  above.  Such  a  nozzle 
produces  a  jet  which  oscillates  from  side  to  side 
and  results  In  large  scale  structures  (shown  In  che 
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Figure  3. 


Schematic  of  the  rotary  valve 
nozzle. 
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Figure  7.  Rotor  jet  flowfield  seen  In  a  moving 
coordinate  ays  ten. 

Another  problem  of  current  Interest  Is  the 
application  of  the  rotor  vortex  generator  to 
reduce  che  drag  of  a  bluff  body  vehicle.  An 
example  is  a  van  in  which  the  objective  i3  to 
reduce  the  drag  but  to  affect  the  volume  of 
the  van  in  a  minimal  way.  Thus,  simply  stream¬ 
lining  the  van  has  severe  limitations.  The 
method  roposed  is  to  employ  the  rotor  to 
enable  che  flow  to  turn  through  a  larger  angle 
and  thus  only  require  a  small  loss  of  van 
volume  due  to  eliminating  che  rear  corner 
as  shown  in  Figure  8.  In  this  figure,  che 
transformation  to  a  coordinate  frame 
structure  dominates  the  flow  and  should  lead 
to  a  drag  reduction  since  it  substantially 
reduces  the  size  of  che  wake  behind  che 
vehicle. 
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chrough  an  analog  to  digical  converter.  The  probe 
positioning  has  been  automated  in  one  direction  and 
is  driven  by  the  microcomputer. 

The  velocity  of  the  air  at  any  given  point  in 
the  flowfield  is  continuously  measured  by  the  hot¬ 
wire  probe.  Since  che  structure  of  che  flowfield 
generated  by  the  rotor  of  Figure  4  is  coherent  and 
synchronously  varying  with  respect  to  che  position 
of  the  rotor  upstream,  che  velocity  is  sampled  with 
respect  to  the  position  of  che  rotor.  Thus  a  com¬ 
plete  "time  picture"  of  che  structure  is  obtained 
as  it  passes  the  probe.  In  addition,  che  coherency 
of  the  flow  structures  allows  che  signals  from 
several  successive  rotations  to  be  averaged  to 
improve  the  signal  co  noise  ratio  and  reduce  che 
effects  of  turbulence  on  the  measurements.  The 
signals  from  the  hot  wire  are  processed  and  recorded 
by  the  microcomputer  through  an  analog  to  digital 
converter. 

The  shaft  of  che  rotor  in  Figure  4  has  been 
fitted  with  a  light  chopper  disk  of  two  cracks. 

One  crack  has  a  single  hole  used  co  indicate  che 
start  of  the  rotor  orbit  and  che  second  crack  has 
thirty-six  equally  spaced  holes  (every  10°)  used  to 
indicate  che  position  of  the  rotor  in  its  orbit. 
These  two  tracks  are  "read"  by  light  emitting  diodes 
(LED),  photo  transistor  pickups.  The  resulting 
signals  are  passed  through  an  interface  which  pro¬ 
vides  signal  conditioning  co  prepare  che  signals 
for  che  microcomputer. 


Figure  3.  Flowfield  produced  by  the  rotor  when 
applied  co  Che  wake  of  a  van. 

A  third  area  of  current  interest  is  che  dy¬ 
namics  of  an  initially  rectilinear  vortex  structure 
near  a  vail  generated  by  che  rotor  device  shown  in 
Figure  4.  A  schematic  of  che  dynamics  which  such 
scruccure  undergoes  is  shown  in  Figure  9.  From  che 
point  of  view  of  application,  such  a  flow  condition 
can  arise  behind  very  rapidly  changed  aerodynamic 
control  surfaces.  The  data  acquisition  system 
described  in  che  following  section  has  been 
developed  specifically  for  chis  investigation  and 
some  cypical  results  will  be  shown  in  chat  section. 
The  results  with  and  without  a  shear  flow  profile 
near  che  wall  are  discussed  in  References  12  and  10 
respectively. 

Instrumentation  and  Data  Acquisition  for  Forced 
Unsteady  Flows 


The  instrumentations  system  has  been  designed 
co  measure  che  instantaneous  velocities  of  c.ne  con- 
vected  coherenc  flow  structures  induced  by  che 
rotor  motion  as  3hovn  in  Figure  4.  The  velocicv 
< two  components)  is  measured  bv  an  "x-wire."  hoc- 
wire  anemometer  and  recorded  by  a  microcomputer 


Figure  9.  Schematic  of  che  deformation  of  an 
initially  rectilinear  vortex. 

The  position  of  che  hot-wire  probe  in  che  flow- 
field  is  partially  compucer  controlled.  The  screatr- 
wise  and  sidewise  positions  of  che  probe  are  manually 
sec  while  che  height  of  the  probe  from  che  tunnel 
floor  is  computer  controlled.  The  probe  is  driven 
up  or  down  by  a  DC  motor  attached  co  a  lead  screw 
wich  a  24  cum-per-inch  pitch.  The  bottom  end  of 
che  lead  screw  is  fitted  wich  a  light  chopper  disk 
with  two  holes.  The  holes  are  read  by  a  LED-phoco- 
cransiscor  pickup.  The  bottom  end  of  che  cravel  of 
che  probe  is  detected  by  another  LED-phototransistor 
pickup  connected  to  a  spring-loaded  probe  placed 
downstream  of  the  hot-wire  probe.  The  signals  from 
the  photo  pickups  are  passed  chrough  the  interface 
electronics  co  the  microcomputer. 

The  hoc-wire  probe  is  positioned  in  che  flow 
ov  first  driving  che  probe  co  boccom  and  then 
bringing  it  up  while  councing  the  turns  of  che  lead 
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The  startup  operation  of  the  unpowered  rotor 
device  is  as  follows:  Even  if  the  rotor  is  well 
balanced  and  relatively  free  to  turn  (i.e.  low 
frictional  resistance),  at  low  freestream  velo¬ 
cities  the  rotor  remains  in  a  static  orienca- 
cion.  As  the  freestreaa  velocity  is  increased,  the 
rotor  assumes  the  position  shown  in  Figure  lib  with 
the  cusp  shape  directly  above  the  shaft.  A  further 
increase  in  the  freestream  velocity  results  in  an 
oscillation  or  rocking  of  che  rotor  about  chis 
position.  The  amplitude  of  che  oscillations  in¬ 
crease  with  freestream  velocity  until  the  rotor 
begins  to  turn  in  the  counterclockwise  direction. 
The  rotor  was  never  observed  to  turn  in  the  clock¬ 
wise  direction. 

sett 


Figure  12a.  Self  driven  rotor  speed  as  a 

function  of  freestream  velocity. 

The  particular  rotor  discussed  in  Reference  10 
has  been  examined  in  detail  to  determine  the 
relationship  becween  che  freescream  velocicy  and 
rotational  speed  it  produces.  The  results,  shown 
in  Figure  12a,  are  well  represented  by  a  linear 
dependence.  This  result  is  emphasized  by  the  non- 
dimensional  results  of  Figure  12b  where  the  rota¬ 
tional  speed  is  reduced  by  V/H;  che  freescream 
velocicy/exposed  rotor  height.  Except  for  low 
freestream  velocities  (where  both  the  velocities 
and  the  rotational  speeds  are  less  accurate),  che 
reduced  rotational  speed  is  nearly  constant. 


Figure  12b.  Self  driven  rocor  speed  non- 

dimensionalized  with  V  4 ;  : ree- 
scream  velocicy/excende-a  r-cor 
height. 


It  should  be  emphasized  that  the  results  of 
Figures  12a  and  b  are  singular  results  for  this 
particular  rotor  shape  and  for  the  frictional 
realstance  of  this  bearing  configuration.  A  change 
in  che  rotor  shape  or  che  frictional  resistance 
would  change  the  value  of  che  curve  but  is  not 
expected  to  affect  the  functional  dependence. 

From  che  point  of  view  of  potential  applica¬ 
tion,  the  fact  that  the  rotor  can  be  employed 
without  the  need  for  an  external  driver  is  signi¬ 
ficant.  It  will  allow  more  widespread  use  of  the 
concept  behind  the  device,  namely  che  potential  use 
of  unsteady  flows  to  achieve  results  superior  to 
chose  found  with  nominally  steady  flows.  All  of  the 
rocor  based  devices  discussed  in  che  preceding 
sections  will  be  positively  affected,  from  che  point 
of  view  of  application,  by  the  possibility  of  self 
induced  rotation. 

Other  Potential  Applications 

A  wide  range  of  potential  applications  exists 
for  unsteady  flow  devices.  Since  an  external  driomr 
is  unnecessary  to  power  the  rocor  vortex  generator, 
the  potential  applications  for  this  particular 
technique  have  improved.  The  following  examples 
are  all  related  to  flow  control  near  a  solid  wall 
and  are  not  meant  co  be  ac  all  inclusive.  Many 
other  embodiments  of  these  concepts  will  appear  in 
che  future. 


For  the  rotor  concept,  a  wall  ;et  version  of 
the  device  is  shown  in  Figures  13a  and  o  where  the 
rotor  has  been  placed  on  opposite  sides  of  a  wail 
;et  exit.  In  either  case,  the  large  scale  vortices 
produced  in  the  jet  will  allow  the  flow  to  remain 
ittacned  to  the  wall  in  situations  where  it  other¬ 
wise  mav  separate.  )f  course,  the  rotor  could  pe 
moiled  on  both  sides  of  the  jet  s imui t aneous Iv  , 
out  it  is  likelv  that  this  will  oe  jnnecessarv  to 
ichieve  tne  lesired  structure. 
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Figure  17.  Orthogonal  descriptors  empLoyed  in 
the  vortex  identification  analysis. 


Figure  IB.  FIowfleLd  aactern  produced  by  :he 
rotacionai.  descriptor. 
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